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Abstract 

Background: Development of durable plant genetic resistance to pathogens through strategies of QTL pyranniding 
and diversification requires in depth knowledge of polygenic resistance within the available germplasm. Polygenic 
partial resistance to Aphanonnyces root rot, caused by Aphanomyces euteiches, one of the most damaging 
pathogens of pea worldwide, was previously dissected in individual mapping populations. However, there are no 
data available regarding the diversity of the resistance QTL across a broader collection of pea germplasm. In this 
study, we performed a meta-analysis of Aphanomyces root rot resistance QTL in the four main sources of resistance 
in pea and compared their genomic localization with genes/QTL controlling morphological or phenological traits 
and with putative candidate genes. 

Results: Meta-analysis, conducted using 244 individual QTL reported previously in three mapping populations 
(Puget X 90-2079, Baccara x PI 180693 and Baccara x 552) and in a fourth mapping population in this study 
(DSP X 90-2131), resulted in the identification of 27 meta-QTL for resistance to A euteiches. Confidence intervals of 
meta-QTL were, on average, reduced four-fold compared to mean confidence intervals of individual QTL. Eleven 
consistent meta-QTL, which highlight seven highly consistent genomic regions, were identified. Few meta-QTL 
specificities were observed among mapping populations, suggesting that sources of resistance are not 
independent. Seven resistance meta-QTL, including six of the highly consistent genomic regions, co-localized with 
six of the meta-QTL identified in this study for earliness and plant height and with three morphological genes 
(At A, R). Alleles contributing to the resistance were often associated with undesirable alleles for dry pea breeding. 
Candidate genes underlying six main meta-QTL regions were identified using colinearity between the pea and 
Medicago truncatulo genomes. 

Conclusions: QTL meta-analysis provided an overview of the moderately low diversity of loci controlling partial 
resistance to A. euteiches in four main sources of resistance in pea. Seven highly consistent genomic regions with 
potential use in marker-assisted-selection were identified. Confidence intervals at several main QTL regions were 
reduced and co-segregation among resistance and morphological/phenological alleles was identified. Further work 
will be required to identify the best combinations of QTL for durably increasing partial resistance to A. euteiches. 

Keywords: Aphanomyces root rot, Pisum sativum, Meta-QTL, QTL diversity. Consistent genomic regions. Resistance 
alleles. Undesirable alleles 
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Background 

Polygenic resistance frequently contributes to partial but 
presumably more durable levels of resistance in culti- 
vated crops [1] and is a major target for many breeding 
programs. Pressure to reduce chemical applications and 
frequent breakdown of major resistance genes in plant 
species encourage the integration of polygenic resistance 
into cultivars of many crops. Polygenic resistance is con- 
trolled by many Quantitative Trait Loci (QTL), which 
express minor to major effects and whose functional 
roles are usually poorly understood [1,2]. The combined 
and diversified use of multiple major or minor-effect re- 
sistance genes should result in both increasing levels of 
resistance [3-6] and in limiting the adaptation of patho- 
gen populations to plant genetic resistances [7-9]. The 
mining of diverse loci and gene functions controlling 
polygenic resistance, combined and deployed together 
with available major resistance genes, is expected to be a 
key strategy for durably improving plant genetic resist- 
ance to pathogens. A depth knowledge of the diversity 
of QTL controlling polygenic resistance among genetic 
resources of cultivated crop species is necessary for 
long-term genetic improvement. St-Clair [10] reviewed 
diversity studies of resistance QTL in several bi-parental 
populations from multiple germplasm sources and the 
subsequent use of marker-assisted-selection (MAS) for 
diverse resistance QTL. Linkage mapping in multi- 
parental population [11] and association mapping in 
connected recombinant inbred lines or natural popula- 
tions [12,13] have been developed for more powerful de- 
tection and comparison of positive alleles among a wider 
diversity of genotypes than in previous studies. Whole 
genome sequencing data in conjunction with genome- 
wide association genetics has been shown a promising 
approach to investigate the diversity of disease resistance 
QTL in plants and to detect with high resolution com- 
mon resistance alleles in a given population [2,14]. QTL 
meta-analysis [15] have also been shown a useful ap- 
proach to obtain a synthetic representation of the diver- 
sity of loci controlling a trait of interest in plants. It 
consists of combining results from previous QTL reports 
into consensus genomic regions, or meta-QTL, associ- 
ated with the trait variation. The position of these meta- 
QTL is estimated with higher accuracy compared to 
position estimates of individual QTL. Statistical methods 
for meta-analysis of QTL have been implemented in two 
main software packages, Biomercator [16] and Meta- 
QTL [17]. This approach has recently been applied in 
plants [18-21] and is especially useful for obtaining an 
overview of loci and for increasing resolution to identify 
candidate genes controlling polygenic resistance [22-29] . 

Aphanomyces root rot, caused by the soilborne 
oomycete Aphanomyces euteiches Drech., is one of the 
most damaging pea diseases worldwide [30] and is a 



major limiting factor to pea production in France. The 
absence of adequate chemical or cultural methods to 
manage the disease requires development of pea cultivars 
with acceptable levels of resistance to sustain pea produc- 
tion. In pea, resistance to A, euteiches has been shown to 
be partial, polygenically inherited [31] and correlated with 
undesirable traits such as long internods, anthocyanin 
production and late-flowering [32] (Roux-Duparque, Pers. 
Comm.). Previous exploration studies of diversity within 
Pisum for resistance to A, euteiches reported that sources 
of resistance were scarce [33-35] (Pilet-Nayel et al unpub- 
lished data). The few sources of resistance identified, were 
integrated into pea breeding programs for root rot resist- 
ance over the past 30 years. These programs resulted 
in the release of germplasm with acceptable agronomic 
characteristics and increased levels of partial resistance in 
the USA [36-42] and in France [43]. In the past 10 years, 
a collaborative French- American research program con- 
ducted genetic dissection of polygenic resistance to A, 
euteiches in four main sources of resistance (90-2079, 90- 
2131, PI180693 and 552) derived from and/or integrated 
into various breeding programs. These sources of resist- 
ance showed the highest and/or most stable levels of re- 
sistance across several different field locations in France 
and in the USA. They presented different genealogies [44], 
except PI180693 and 90-2131 which were related [38] but 
differed in their levels of resistance. Results from this col- 
laborative program have been reported for three mapping 
populations [31,45,46]. From a RIL population derived 
from the cross Puget (susceptible to A, euteiches) x 90- 
2079 (resistant to A, euteiches in the USA), 14 QTL were 
associated with Aphanomyces root rot resistance, includ- 
ing three QTL {Aphl, Aph2 and Aph3) that were detected 
over various environments and isolates [31,45]. In two 
RIL populations derived from crosses between the suscep- 
tible variety Baccara and the resistant germplasm lines 
PI 180693 and 552, a total of 23 genomic regions control- 
ling partial resistance to A. euteiches were identified. 
These regions included five highly stable regions across 
environments (France and USA), isolates and genetic 
background, as well as 13 epistatic interactions [46]. QTL 
results obtained from a fourth RIL population derived 
from the cross between DSP (susceptible) and 90-2131 
(resistant) are reported in this manuscript along with QTL 
comparisons between the independent mapping studies. 
This comparison will provide a more complete view of the 
diversity of Aphanomyces resistance loci available from 
the four sources of resistance. Greater understanding of 
the resistance QTL will allow QTL pyramiding and diver- 
sification strategies to be developed, in an effort to prevent 
the pathogen from overcoming the crucial and limited 
level of resistance. Comparisons between Aphanomyces 
resistance QTL and QTL controlling other agronomic 
traits will aid selection for positive associations and avoid 
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selection of negative linkages between resistance and un- 
desirable alleles. 

The objectives of this study are (i) to provide an 
overview of the diversity and the structural genomic 
organization of QTL involved in resistance to A. euteiches 
in pea, by compiling QTL mapping data obtained from 
four main sources of resistance, (ii) to identify positive 
and negative co-segregations between Aphanomyces re- 
sistance alleles and alleles controlling morphological and 
phenological traits and (iii) to identify putative candidate 
genes underlying reduced confidence intervals of main 
meta-QTL using pea-M truncatula translational genom- 
ics. We conducted a meta-analysis of QTL for resistance to 
A, euteiches in pea integrating (i) QTL previously detected 
for resistance to A, euteiches from three resistance sources 
(90-2079, 552, PI180693) [31,45,46] and (ii) QTL reported 
for the first time in this study from the additional source of 
resistance, 90-2131. We also performed a meta-analysis of 
QTL controlling morphological and phenological traits 
(plant height and earliness), from phenotypic data and 
QTL results obtained in three of the four RIL populations 
studied for Aphanomyces resistance. Finally, we compared 
genomic localizations of the identified morphological and 
phenological meta-QTL to the ones of Aphanomyces 
resistance meta-QTL. We identified candidate genes 
underlying main resistance meta-QTL using the pea- 
M. truncatula translational toolkit of Bordat et al, [47]. 

Results 

Genetic analysis of resistance to A, euteiches in the DSP x 
90-2131 RIL population. 

The genetic map, constructed from the DSP x 90-2131 
RIL population, was comprised of a total of 168 markers, 
including 107 SSRs, 56 RAPDs, three genes of known- 
function (two flowering related genes, LD and FPA [48], 
and a putative sugar transporter, SugTrans [49]) and two 
morphological markers (r and PI) distributed over nine 
linkage groups (LG) (Figure 1). The map covered 1046 cM 
Kosambi, corresponding to about 77% of the pea SSR 
reference genetic map [50]. Five of the 168 markers (3%) 
showed significant distorsion from the expected Mendelian 
ratio (a=0.01). Of the 168 markers, 86 (51%), equally 
distributed over nine LG, were common to the pea refer- 
ence genetic map [50], including 74 SSR, 11 RAPD and 
one morphological marker {PI), These common markers 
were ordered with high colinearity between the two genetic 
maps, except in the distal part of LGIIb, where the order of 
a block of markers was inverted. Compared to the refer- 
ence map, the distal part of LGI and the pericentromeric 
regions of LGII and LGIII remain uncovered by markers 
on the DSP x 90-2131 map. 

Analysis of the DSP x 90-2131 RIL resistance data in 
controlled conditions showed that the five pea lines used 
as controls ranked as expected for root rot index (RRI) 



scores with all but one of the six A, euteiches strains, 
confirming the strain patho types as described in [51] 
and [52] (Additional file 1). RRI scores in the DSP x 90- 
2131 RILs showed highly significant genotypic and block 
effects {P < 0.001) for each strain with the exception 
of Ae85. Mean-based heritability of RRI was moderate 
(/z^ = 0.55 with Ae85 strain) to high (/z^ = 0.91 with RB84 
strain). Distribution of RRI adjusted means in the RIL 
population tended to fit a normal curve with each of the 
six strains except RB84 (Additional file 2). Lower levels 
of symptoms were observed in the RIL population with 
the Ae85 strain than with the five other strains {[i = l,9 
and 2.9-3.6, respectively). Transgressive RIL segregants 
with increased resistance and susceptibility compared to 
parental values were observed with each of the six 
strains. 

Analysis of the DSP x 90-2131 RIL resistance data in 
each of the 11 field environments studied showed highly 
significant genotypic effect (P<0.002) for all but two vari- 
ables (RI-2002-RRI, RI-2003-ADI1) and significant block 
effects (P<0.002) for half of the variables. Mean-based 
heritabilities of traits, ranging from 0.25 (RI-2003-ADI1) 
to 0.87 (PLM-2000-ADI1), were low {h^ < 0.40) for five 
variables, moderate (0.40 <h^ < 0.70) for nine variables 
and high (h^ > 0.70) for six aerial decline index (ADI) 
variables (Additional file 3). Frequency distribution of 
the adjusted mean scores for each variable in the RIL 
population tended to fit normal curves, with the excep- 
tion of the ADI scores obtained in Pullman in 2000 and 
2003 (Additional file 3). Differences were observed be- 
tween means and ranges of the RIL population for RRI 
and ADI variables, depending on the location and year. 
In all experiments, transgressive segregants, either more 
susceptible or more resistant than the parents, were 
observed. 

Half of the DSP x 90-2131 RIL phenotypic data from 
different field criteria and environments were signifi- 
cantly {P < 0.001) and positively correlated with most of 
the other field data (Additional file 4). Data with low 
heritability values were mostly poorly correlated with the 
other field data (RI-2002-RRI, RI-2003-ADI1, DI-2003 
-RRI and LS-2003-ADI1). Data from evaluations in con- 
trolled conditions for the six A, euteiches strains were 
highly correlated (P < 0.001) with each other. The RB84 
strain data was also positively and significantly corre- 
lated with most of the RRI and ADI scores assessed in 
French and American field conditions, contrary to the 
other strain data which were significantly correlated with 
few or no field scoring data. 

A total of 79 additive-effect QTL, corresponding to 25 
genomic regions distributed over nine linkage groups, 
were detected in the DSP x 90-2131 RIL population for 
Aphanomyces root rot resistance evaluated in 11 field en- 
vironments and with six strains in controlled conditions 
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(See figure on previous page.) 

Figure 1 Genomic localization of QTL for resistance to A. euteiches in the DSP x 90-2131 RIL population. Genetic map constructed from 
1 1 1 RILs derived from tlie cross DSP x 90-2131 and genomic localization of additive-effect QTL detected for Aphanomyces root rot resistance in 
1 1 field environments over four years and five locations and against six strains (RB84, Ael06, Ae85, Ae87, Ae78 and Ael09) in controlled 
conditions, using two resistance criteria (Root Rot and Aerial Decline Indexes, RRI and ADI, respectively). Linkage groups (LG) assigned to the SSR 
pea reference map (Loridon et oi, 2005) are named from I to VII. The size of each LG is indicated in cM Kosambi. Marker names are indicated on 
the left of each LG. Markers with biased allelic segregation are indicated by one asterisk. Markers common to the SSR pea reference map (Loridon 
et oi, 2005) are underlined. Lengths of additive-effect QTL boxes correspond to the LOD-1 support interval from the peak marker. 



(Figure 1, Additional file 5). A major QTL, AePs7.6, previ- 
ously identified in [46], was consistently detected on 
LGVII from ADI and/or RRI data collected in the 11 field 
environments studied and for five of the six strains in con- 
trolled conditions. The Ae-Ps7.6 genomic region covered 
about 35 cM and accounted for the majority of the pheno- 
typic variation, yielding 59.8% for resistance to the RB84 
strain. Two QTL, Ae-Ps3,l and Ae-PsSJ, previously iden- 
tified [46], were also consistently detected on LGIIIa and 
LGV (near to the R locus), especially from ADI and/or 
RRI field data (Ae-PsSJ) and from the six strain data in 
controlled conditions (Ae-PsS.l). Five genomic regions on 
LGIV and LGVI {Ae-Ps4.1, Ae-Ps6.h Ae-Ps43, Ae-Ps4A 
and Ae-Ps6,4), were consistently associated with three to 
five individual field and/or controlled condition resistance 
QTL. Seventeen QTL, distributed on the seven LG, were 
less consistently identified from one or two variable (s) of 
the same criterion. For the eight total genomic regions 
detected from at least three variables, no consistent 
specificities were observed across either environments or 
strains, except for the three regions that were more 
specifically identified from French field {Ae-Ps4.1 and Ae- 
Ps6,l) or controlled condition {Ae-PsS,l) data. Four of the 
eight regions {Ae-Ps7,6, Ae-PsS.l, Ae-Ps4.4 and Ae'Ps4,3) 
explained more than 15% of the phenotypic variation 
observed for at least one variable. Resistance was contrib- 
uted by 90-2131 and DSP alleles, whatever the variable, at 
four {Ae-Ps7.6, Ae-PsSJ, Ae-Ps3.h Ae-Ps4J) and two 
{Ae-Ps4.4, Ae-Ps6.1) of the eight regions, respectively. Four 
significant pairwise epistatic interactions were identified 
for increasing resistance to A. euteiches, especially to the 
RB84 strain. All the interactions involved at least one 
minor QTL, especially A 7. 3 (Additional file 6). 

Genetic analysis of earliness and plant height from DSP x 
90-2131, Baccara x PI1 80693 and Baccara x 552 RIL 
populations 

Analysis of variance of earliness and plant height data in 
each environment and population showed a significant 
genotypic effect (P< 0.0001) and no significant block 
effect {P > 0.05) for all variables. Mean-based heritabilities 
of the flowering traits (FLO) and plant height variables 
(HT) evaluated in a three-block design were very high 
(/z^>0.80), except for FLOl assessed at Dijon in 2003 
(Additional file 7). Frequency distributions of the FLO and 



HT scores closely followed normal curves in the DSP x 
90-2131 and Baccara x 552 populations. FLO distribu- 
tions tended to fit bimodal curves in the Baccara x 
PI 180693 population in most of the environments tested, 
despite the similar values of parental lines. The FLO and 
HT scores were significantly (P< 0.001) and negatively 
correlated with no, few and most Aphanomyces resistance 
variable scorings in the Baccara x 552, DSP x 90-2131 
and Baccara x PI 180693 populations, respectively (data 
not shown). FLO and HT values were highly correlated in 
each population, across all environments (infested/healthy). 

A total of 33 additive-effect QTL corresponding to 14 
genomic regions were detected for the flowering traits 
(FLOl, FL02), from the three RIL populations studied. 
Additionally, three additive-effect QTL corresponding to 
three genomic regions were detected for plant height 
(HT) (Additional files 8 and 9, and Figure 2). No co- 
localization was observed between genomic regions 
detected for FLO and HT traits. Four genomic regions 
were consistently detected for FLO, each containing four 
to six individual QTL. Three of them {Flo-Ps2.2, Flo-Ps3.1, 
FlO'Ps72) explained up to 40.8% of the phenotypic vari- 
ation while FlO'Psl.2 explained up to 10% of the variation. 
Two other regions (Flo-Psl.l and Flo-Ps6.3) were also 
identified, each containing two individual QTL. Flo-Ps6.3 
explained up to 40% of the phenotypic variation. Among 
the six regions, alleles for earliness at Flo-Psl,l, Flo-Ps2.2 
and Flo-Ps3,l were contributed by Baccara, while they 
were contributed by 90-2131 at Flo-Ps6.3 and PI180693 
or 552 at Flo-PsL2 and Flo-Ps7.2. Eight FLO QTL with 
minor effects {R^ = 3.5 to 16.3%), were specifically detected 
in one environment. The three QTL identified for plant 
height were detected with minor-effects, and at these loci, 
the 90-2131 alleles contributed to shorter plants. 

QTL meta-analysis and co-localisations between 
Aphanomyces resistance and morphological/phenological 
meta-QTL 

QTL meta-analysis of Aphanomyces resistance and mor- 
phological/phenological traits was performed using three 
mapping studies in the four pea RIL populations de- 
scribed in Table 1. The QTL mapping study conducted 
using the Puget x 90-2079 RIL population [31,45] was 
updated by adding 53 markers (41 SSRs and 12 RAPDs) 
to the genetic map. Thirty five of the added markers were 
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Table 1 QTL mapping populations and studies used for meta-analysis 



Cross Populations size. Donor Genetic map Phenotypic evaluation QTL detection Reference 

type and generation method^ 



ruget X yu-zU/y 


1 T7 DM /CI n^ 
1 z/ KIL (r 1 Uj 


1 ICPi A A DC 

(USA) 


324 markers 


Resistance to A. euteiches 

Roci c"l"3 nr"o fr\ A oiifoirhtOQ 
ncblbLdI ILc LU A\. cULclLilcb 


LIIVI 


riiet-Nayei 
etal. [31,45] 


Baccara x 
PI 180693 


178 RIL(F8) 


INRA (France) 


224 markers 
(consensus genetic 
map) 


Morpliological and 
plienological traits 

Resistance to A. euteiches 


CIM 


Hamon 
et al. [46] 


Baccara x 552 


178 RIL(F9) 


INRA (France) 




Morpliological and 
plienological traits 
Resistance to A. euteiches 

Resistance to A. euteiches 
resistance 


CIM 


Hamon 
et al. [46] 


DSP X 90-2131 


111 RIL(FIO) 


USDA-ARS 
(USA) 


168 markers 


Morpliological and 
plienological traits 


CIM 


This study 



^ Composite interval mapping. 



common to at least one of the two other individual genetic 
maps used The updated Puget x 90-2079 genetic map, 
covering 1523 cM, comprised a total of 377 markers (data 
not shown). Individual QTL detected on the previous map 
[31,45] were confirmed and the updated QTL information 
were used for the meta-analysis. 

Consensus map 

The consensus marker map was established from the 
three individual genetic maps (Puget x 90-2079, DSP x 
90-2131, Baccara x PI180693/552) using the Meta-QTL 
software. The map comprised a total of 619 markers 
(31% SSR; 33% RAPD; 32% AFLP; 2.4% STS, 1.6% iso- 
zyme or morphological markers) and covered 1513 cM 
Kosambi (Additional file 10). Eighty-six percent of the 
markers mapped to only one of the individual maps 
studied, while 14% mapped to two or three individual 
maps with consistent marker order and regular distribu- 
tion on all LG (8 to 22 common markers per LG) with 
exception of LGI. Fifteen percent of the 619 markers of 
the consensus map were also mapped on the reference 
genetic map described by [53] and had consistent posi- 
tions between the two maps. 

Meta-QTL for resistance to A. euteiches 

A total of 244 individual additive-effect QTL detected for 
resistance to A. euteiches in the three mapping studies 
were projected onto the consensus marker map (Figure 2). 
Projection of all individual QTL allowed three to five main 
genomic regions per LG to be identified, each consisting 
of individual QTL with overlapping confidence intervals. 
Removal of redundant or single-individual QTL left 115 
independent QTL for the meta-analysis. Meta-analysis 
identified 27 meta-QTL for resistance to A, euteiches dis- 
tributed over the seven LG, with three to four meta-QTL 
per LG (Table 2, Figure 2). Each meta-QTL was 



comprised of 1.2 to 11.8 independent QTL (mean = 4.2 ± 
2.5 independent QTL per meta-QTL) corresponding to 
one to 41 initial individual QTL (mean = 9 ± 8 initial QTL 
per meta-QTL). Confidence intervals of meta-QTL ranged 
from 0.9 to 53.1 cM (Kosambi), with a mean of 13.7 ± 
12 cM, and were lower than 20 cM for 85% of the meta- 
QTL. The average width of the independent QTL confi- 
dence intervals was reduced from 0.9-fold to 28.4-fold 
(4.1 ± 5.2 in average), depending on the meta-QTL. Among 
the 27 meta-QTL identified, 11 were highly consistent, 
including eight meta-QTL comprised of at least 12 initial 
individual QTL and three meta-QTL comprised of eight to 
11 initial QTL (Table 2). The remaining 16 less consistent 
meta-QTL aggregated individually between one and six 
initial QTL. The 11 consistent meta-QTL corresponded to 
ten independent genomic regions. In seven of these re- 
gions, initial QTL were clustered in highly consistent posi- 
tions, namely on LGI {MQTL-Ae3\ LGII {MQTL-AeS\ 
LGIII {MQTL-Ae8/Ae9\ LGIV {MQTL-Ael2 and MQTL- 
AelS\ LGV {MQTL'Ael6/Ael7) and LGVII {MQTL'Ae2S/ 
Ae26), Eight of the 11 consistent meta-QTL aggregated ini- 
tial QTL detected from either three or four mapping popu- 
lations, and all but one of the remaining meta-QTL were 
detected from two mapping populations. At each meta- 
QTL that contained initial QTL from at least two mapping 
populations, two to five parental alleles, mostly from the 
resistant parent, contributed to resistance with a mean of 
3.2 alleles per meta-QTL. Five of the 24 meta-QTL com- 
prised of initial QTL from at least two field environments 
were specific to either French {MQTL-Ael2, MQTL-Ae20 
/Ae21) or USA {MQTL-AelS, MQTL-AelT) environments, 
while the remaining 19 were not environment-specific. All 
meta-QTL including initial QTL detected from controlled 
condition data co-located with initial QTL identified from 
at least one field environment. None of the meta-QTL 
comprised of a minimum of two initial controlled 



Table 2 Meta-QTL detected for resistance to A. euteiches and morphological/phenological traits 
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Table 2 Meta-QTL detected for resistance to A. euteiches and morphological/phenological traits (Continued) 



Morphological/Phenological traits 



I MOTL- 
Morphol ' 

MQTL- 
Morpho2 ' 

II MOTL- 
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^ Meta-QTL names are coded as follows : MQI/.-zAe = Meta-QTL detected for resistance to A euteiches; MQTL-Morpho - Meta-QJl detected for morphological and phenological traits (flowering time and plant height). On 
each Linkage Group (LG), MQTL-Ae and MQTL-Morpho, which co-localize, are mentioned by asterisks. Meta-QTL indicated in bold correspond to the seven main consistent genomic regions identified for resistance to 
A. euteiches. 

^ QTL position at the LOD (Log of likelihood ratio) score peak, from the first marker of the linkage group in cM Kosambi. 
^ 95% confidence interval of the meta-QTL given by MetaQTL software, in cM Kosambi. 

^ Mean and standard deviation of confidence intervals (estimates according to [55]) of all independent QTL belonging to each meta-QTL, in cM Kosambi. 
^ Number of initial QTL (redundant and independent) underlying each meta-QTL. 

^ Range of percentage of phenotypic variance (R^) explained by individual initial QTL underlying each meta-QTL. 
^ Range of LOD-1 confidence intervals of all initial QTL underlying each meta-QTL. 

^ Number of variables for which each RIL parental allele contributed to increased resistance to A. euteiches, increased earliness or short plants, at each meta-QTL. No values indicate missing data. 
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condition QTL, were specific to a single strain. Eighty per- 
cent of the meta-QTL were not specific to root or aerial 
scoring criteria used in the different studies, even if some 
of them included initial QTL that were primarily detected 
from root criteria {MQTL-Ael6) or aerial criteria (MQTL- 
Ae27). 

Meta-QTL for morphological/phenological traits and co- 
localizations with Aphanomyces resistance meta-QTL 

Projection of 36 initial additive-effect QTL detected for 
earliness and plant height onto the consensus map outlined 
15 distinct genomic regions with QTL clusters, including 
21 redundant or single initial QTL. Meta-analysis of the 
remaining 15 non-redundant initial QTL resulted in the 
identification of seven meta-QTL distributed over six LG 
(Table 2, Figure 2). Each meta-QTL was comprised of one 
to three independent QTL (mean = 2.0 ±0.7 independent 
QTL per meta-QTL), corresponding to two to six initial 
individual QTL. Confidence intervals of the meta-QTL 
ranged from 0.7 to 28 cM and were reduced by 1 to 37.1 
times compared to mean confidence intervals of clustered 
initial QTL. The seven meta-QTL corresponded to the six 
most consistent FLO genomic regions {Flo-Psl.l, FloPsl.2, 
Flo-Ps2J/2.2, Flo-Ps3.l Flo-Ps6.2/6.3, Flo-PsZ2), and to a 
cluster of two FLO and one HT overlapping QTL {MQTL- 
MorphoS), Five of the seven meta-QTL were identified 
from at least two populations and from diverse en- 
vironments. 

Based on the 33 meta-QTL that contained more than 
one initial QTL for either Aphanomyces resistance or for 
plant earliness and height, six genomic regions were 
detected with overlapping intervals of resistance and 
phenological QTL (Figure 2, Table 2). In the regions of 
meta-QTL co-localization, nine negative and five positive 
associations were found between parental alleles contrib- 
uting to resistance and early flowering and/or short stat- 
ure. Positive associations involved 552, PIl 80693 and DSP 
alleles. PIl 80693 and/or 552 alleles contributing to resist- 
ance and normal leaf type or anthocyanin production were 
in coupling phase in meta-QTL regions on LGI (MQTL- 
AeS) and LGII {MQTL-AeS/Ae6). 90-2079 and 90-2131 
alleles contributing to resistance and round seed type were 
also in coupling phase in the meta-QTL region on LGV 
{MQTL-Ael6/Ael7), 

Identification of candidate genes underlying meta-QTL 

Using the V^di-Medicago truncatula translational toolkit 
of Bordat et al, [47], a Ust of genes among 5460 pea 
Unigenes was identified as positional candidates within 
the support interval of all but one {MQTL-AelS) of the 
seven highly consistent resistance meta-QTL regions. 
Adjacent markers to the meta-QTL, which were com- 
mon between the consensus marker map developed in 
this study and the pea functional map of [47], were used 



to define intervals containing candidate genes 
estabUshed in [47]. A total of 318 genes underlying main 
meta-QTL were identified, including 14 to 91 genes per 
meta-QTL (Additional file 11). Out of them, 264 were 
similar to genes with a putative function, among which 
17% had a known-function related to plant disease re- 
sistance. These genes corresponded to resistance gene 
analogs (RGA) or genes involved in plant-pathogen rec- 
ognition (LRR, Leucine zipper proteins), genes involved 
in defense or signal transduction (protein kinases, en- 
zymes of oxidative stress and cellulose synthesis, 
pathogenesis-related-proteins, transcription factors, 
heat-shock proteins, cell division proteins, cyclin-like F- 
Box proteins) and genes involved in resistance to other 
diseases. Particularly, in meta-QTL regions with smallest 
confidence intervals, genes coding for LRR proteins 
{MQTL-Ae3) and protein kinases {MQTL-Ae8/Ae9, 
MQTL-Ael2) could be suggested as good candidates. In 
the other meta-QTL regions, clusters of genes coding 
for heat shock proteins {MQTL-Ael6/Ael7) and a plant 
disease response protein-encoding gene {MQTL'Ae26) 
were also striking putative candidates to the resistance. 

Discussion 

A moderately low diversity of loci controls quantitative 
partial resistance to A. euteiches in four main sources of 
resistance in pea 

This study is the first that provides an overview and 
comparison from multiple studies of genetic loci con- 
trolling resistance to A. euteiches in pea. It is one of the 
first reports developing a comprehensive picture of the 
genetic architecture of disease resistance in cultivated le- 
gumes [24,54]. Similar studies have proven valuable in 
other plant species for determining key resistance loci 
useful in MAS and identifying candidate genes under- 
lying resistance QTL [22,23,26]. 

In this study, 244 individual QTL detected for 
Aphanomyces resistance from eight variables collected 
in four pea RIL populations over a total of 29 field envi- 
ronments (two-five USA-FR locations over 10 years) 
and 12 controlled condition assays (two-six strains per 
population) were analyzed for meta-analysis. A total of 
27 meta-QTL were identified, including 11 consistent 
meta-QTL each containing more than eight initial indi- 
vidual QTL. Seven main genomic regions clustering nu- 
merous initial QTL with high position consistency were 
highlighted. The number of consistent meta-QTL iden- 
tified is quite low given the four diverse partially resist- 
ant pea lines studied and the complex inheritance of 
partial resistance to A, euteiches. 

The meta-QTL analysis conducted in this study accur- 
ately compared genomic positions of individual QTL iden- 
tified from different studies and refined the confidence 
intervals of the main genomic regions associated with 
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resistance. The consensus map preserved marker order of 
individual maps based on 14% of common markers, and 
of the reference genetic map [53] based on 15% of com- 
mon markers. Identification of over-reduced meta-QTL 
confidence intervals was minimized by clustering individ- 
ual QTL using non-redundant independent QTL with 
confidence intervals estimated from their values [55]. 
This approach resulted in the identification of three to 
four meta-QTL per LG, and their confidence intervals 
were reduced four-fold, on average, compared to mean 
confidence intervals of the initial independent QTL. Con- 
fidence interval reduction was correlated to the number of 
clustered independent QTL available for refining each 
meta-QTL (r^ = 0.21), as shown by the smallest confidence 
intervals (0.9-6.9 cM) observed for the eight most consist- 
ent meta-QTL. 

In contrast to previous similar meta-QTL studies of 
quantitative resistance [25,26], few meta-QTL specificities 
were observed associated with the parental source of re- 
sistance. All but one meta-QTL clustering at least two of 
the initial QTL were detected in two, three or four map- 
ping populations. Each of these meta-QTL comprised of 
two to five alleles that contributed to the resistance. At a 
given meta-QTL, meta-analysis did not provide effect esti- 
mates of the different alleles, i.e. classification of alleles 
according to their effect levels. However, at several meta- 
QTL, variations in additive effects and were observed 
between initial QTL detected from different RIL popula- 
tions. Initial QTL detected from the DSP x 90-2131 
population at MQTL'Ae26 had higher average additive ef- 
fects and R^ values (a = 0.3, R^ = 20%) than those detected 
from the Baccara x PI180693 (a = 0.21, R^ = 10%) and 
Baccara x 552 (a = 0.14, = 10%) populations. This sug- 
gests that some parental alleles contributing to the resist- 
ance, generally detected with high consistency, could have 
stronger effects than others. Multiple parental alleles asso- 
ciated with resistance at one meta-QTL could correspond 
to identical or different alleles of a single gene or to clus- 
ters of closely linked genes. Despite the reduced confi- 
dence intervals obtained for several meta-QTL, fine 
mapping of meta-QTL will be required to determine the 
correct allele/gene hypothesis. According to pedigree in- 
formation available in the literature, the four resistant 
germplasm studied are not independent and may share 
common resistance genes, which is supported by the low 
specificity of most of the meta-QTL identified in the RIL 
populations. PI 180693 was one of the parents used in 
crosses to develop 90-2131 but was not reported in the 
genealogy of 90-2079, which derived its resistance from 
MN313 [36,38]. Little information is available about the 
parentage of 552, which was derived from several recur- 
rent selection cycles conducted by Lewis and Gritton [56] . 
The literature contains no data about genetic distances 
between the sources of resistance to A. euteiches used. 



However, in our study, 20 meta-QTL underlied different 
combinations of resistance alleles from two of the four re- 
sistant parent, frequently including 552 alleles, and nine 
meta-QTL included resistance alleles derived from three 
of the four sources of resistance. These results suggest that 
a common genetic background was used in the different 
Aphanomyces resistance breeding programs which se- 
lected the different sources of resistance studied. 

Specificity of Aphanomyces resistance meta-QTL to en- 
vironments (years, locations in France and USA), strains 
and scoring criteria, was not observed, as discussed in [46] 
for QTL identified in Baccara x PI 180693 and Baccara X 
552 populations. In this study, we confirmed the low spe- 
cificity of resistance QTL in the DSP x 90-2131 popula- 
tion, especially at the two highly consistent meta-QTL 
MQTL-Ael6 and MQTL-Ae2S/Ae26, towards a number 
of French and USA environments or a number of A. 
euteiches strains from the same or different pathotypes 
tested than in [46]. However, some genomic regions were 
still observed to be highly specific to French environments 
{MQTL-Ael2, MQTL-AeS), controlled condition tests with 
different strains {MQTL-Ael6) or aerial scoring criteria 
{MQTL-Ae27). 

Aphanomyces resistance alleles at most the consistent 
QTL co-segregate with alleles at genes/QTL controlling 
morphological and phenological traits 

In this study, we identified co-segregating alleles contrib- 
uting to resistance and morphological or phenological 
traits in seven genomic regions, including all but one 
(MQTL-Ael2) of the seven main consistent genomic re- 
gions associated with resistance to A. euteiches and one 
less consistent resistance meta-QTL (MQTL-Ael). In 
these genomic regions, alleles derived from resistant par- 
ents were often associated with undesirable alleles for 
dry pea breeding (late-flowering/higher plant, normal 
leaves (4/), colored flowers {A); Table 2). Co-segregation 
between resistance and favorable alleles (round seeds 
{R), early-flowering plant), especially ones derived from 
552, were identified in some genomic regions {MQTL- 
Ae3, MQTL-AelS, MQTL-Ael6/Ael7). Marx [32] also 
showed that tolerance to A. euteiches was genetically 
associated with several wild- type alleles, including A 
(colored flowers) which was found to be associated with 
MQTL-AeS in this study {Ae-Ps2.2 in [46]), as well as Le 
(tall plants) and PI (black hilum of the seeds). Le and PI 
were not associated with resistance in this study, 
according to results obtained from one of the RIL po- 
pulations studied (Baccara x PI180693, [46]). We previ- 
ously reported that R (round seeds) and Af (normal 
leaves) were other alleles associated with Aphanomyces 
resistance at QTL Aph2 and Ae-Psl.2, i.e. AphS, in pea 
[31,46]. In this study, we demonstrated again the positive 
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co-segregation between R and resistance alleles in the 
DSP X 90-2131 RIL population at the Ae-Ps5J QTL, i.e. 
Aph2, Fondevilla et al [57] and Prioul et al [58] previously 
reported QTL for earliness {dflLhflol and dfiILZflo2) that 
we localized in same genomic regions as two main consis- 
tent QTL identified in this study for earliness and resist- 
ance to A. euteiches {Flo-Ps2,2/MQTL-AeS on LGII and 
Flo-Ps3J/MQTL-Ae8/9 on LGIII). Co-segregation between 
alleles associated with late-flowering and partial resistance 
to Mycosphaerella pinodes were identified in the region on 
LGIII [54,55] and the same associations were observed in 
this study for resistance to A, euteiches. 

Diverse candidate genes underlie Aphanomyces 
resistance meta-QTL 

The negative or positive associations observed in this 
study between resistance and morphological/phenological 
alleles may correspond to pleiotropic genes controlling 
plant architectural or developmental traits or to different 
closely linked genes. Pleiotropic genes controlling plant 
architecture or development have already been suggested 
as good candidates for underlying resistance QTL [2]. In 
pea, clusters of QTL controlling numerous traits related 
to plant morphology, seed protein content and yield, ni- 
trogen nutrition, and root architecture were mapped close 
to architecture and developmental genes, especially Le 
and Af, which were shown to be localized in genomic 
regions having pleiotropic effects [59,60]. Our study sup- 
ports the hypothesis that A/ is a good candidate gene for 
pleiotropy, impacting Aphanomyces resistance (at MQTL- 
AeS) and earliness. Evaluation of NILs or mutants at the 
A/gene [60] for resistance to A, euteiches will be useful for 
validating the pleiotropy hypothesis. 

Clusters of different closely linked genes, especially re- 
sistance genes [61], located in repetitive non-coding se- 
quences have also been reported in plants, with increasing 
plant genome sequences available [62,63]. In our study, 
several Aphanomyces resistance QTL were localized in 
regions previously reported to contain QTL controlling 
resistance to other stresses. Based on comparative map- 
ping of SSR markers, the MQTL-Ae8-9 {Ae-Ps3J) region 
coincides with a main consistent resistance QTL to M. 
pinodes {MpIILS, [57,58]), a major frost tolerance QTL 
co-localizing with the Hr (photoperiod high-responsive 
flowering) locus [64] and a minor QTL for resistance to 
Fusarium oxysporum race 2 {FwnS.l, [65]). Similarly, the 
MQTL-AelS (Aphl), MQTL-Ae2S/Ae26 {Ae-Ps7.6a) and 
MQTL-Ae20 {Ae-Ps6.1) were identified in similar regions 
as the major QTL Fwn4,l for resistance to F, oxysporum 
race 2 [65], a major resistance QTL to Fusarium solani 
[66], along with the Qrw/' QTL for resistance to Uromyces 
fabae [67] and the QTL FRR3 for resistance to F, 
solani [68], respectively. Regions of resistance QTL co- 
localizations have been suggested to underlie clusters of 



resistance gene analogs [2]. Especially, the main-effect 
MQTL-Ae2S/Ae26 region has been located in a region of 
RGA clusters [69], suggesting that defeated resistance 
gene(s) may underlie the QTL. However, other mecha- 
nisms may also be suggested for resistance QTL, i.e. 
defense, signal transduction, as highlighted by the diversity 
of candidate genes identified in this study for main meta- 
QTL. Our results do not favor one particular hypothesis 
for molecular basis of resistance QTL rather than another, 
corroborating Ballini et al and Danan et als conclusions 
[22,23]. Fine mapping and mutants studies, as well as 
genome sequencing efforts, wiU be necessary to discover 
causal genes underlying resistance meta-QTL to A. 
euteiches in pea and to validate hypotheses currently pro- 
posed in this study regarding co-localizations between 
candidate genes and resistance QTL. 

Conclusion 

This study describes alleles that significantly contribute 
to resistance to A. euteiches and their positive or nega- 
tive associations with morphological/phenological traits. 
From the meta-analysis conducted in this study, a choice 
of alleles at meta-QTL corresponding to seven highly 
consistent genomic regions controlling resistance can be 
recommended, with their associated markers, for use in 
MAS (Table 3). Resistance alleles at the most consistent 
region, MQTL-Ae25/Ae26 on LGVII, and especially the 
high-effect allele from 90-2131, appear to be the best 
choices for improving resistance. Resistance allele de- 
rived from 90-2079 at MQTL-AelS {Aphh [31]) will be 
useful only in some USA environments from which the 
QTL was specifically detected. Resistance aUeles at the 
five other genomic regions could be prioritized for use 
in MAS, depending on selection focus and breeding 
objectives {i.e. multiple-environment effects, markers, 
confidence intervals, association with unfavorable al- 
leles). RILs from the DSP x 90-2131 population carrying 
90-2131 alleles at most of markers linked to the main 
meta-QTL on LGIII, LGV and LGVII (Table 3) have 
been recently released [42] and wfll be useful in breeding 
for resistance to A. euteiches. 

Because sources of resistance are limited [35] and the 
number of highly consistent resistance QTL highlighted 
from this study is moderately low, the most consistent 
genomic regions controlling resistance should be carefully 
managed to prevent reduced effectiveness due to potential 
development of new virulent isolates of A. euteiches. In al- 
falfa, strains of A. euteiches from race 2 have been 
reported to overcome resistance to race 1 in cultivars 
grown in the USA [70] . In Medicago truncatula, the level 
and nature of the effect of a major Aphanomyces resist- 
ance QTL {AERl/prAel, [71,72]) have been shown to vary 
depending on the pathotypes of A. euteiches strains [73], 
suggesting that QTL might have variable effects on 



Table 3 Information about the seven highly consistent genomic regions useful for MAS 



Meta-QTL 
(other names) 


LG 


Consistent 
resistance 
enhancing allele 


Number of variables 
for which allele 
contributed to 
resistance 


Effects on 

environments, strains^ 


min-max^ 


Confidence 
interval (cM)*^ 


Recommended markers 
for MAS associated with 
resistance allele*^ 


Unfavorable 

morphological/ 

phenological 

alleles co-segregating 

with resistance alleles 

detected in this study 


MQTL-Ae3 {Ae-Psl .2, AphS) 


1 


552 

PI 180693 
90-2079 


7 
5 
3 


Field (RI-DI), RRI-ADI 
Field (RI-DI-LS), RRI-ADI 
CC (SP7, Ael06) 


8- 14% 
4-12% 

9- 15% 


4-55 


D2] (270), AFO 16458 
(180), AC75 (280) 


normal leaves 
normal leaves 
normal leaves 


MQTL-Ae5/Ae6 {Ae-Ps2.2) 


II 


PI 180693 
552 


8 
3 


Field (RI-DI), ADI-DW CC 
(Ael09) 

Field (Dl), RRI CC (RB84, 
Ael09) 


7-27% 
6-8% 


35-44 


AA372.1 (260;, AB112 (390), 
AD83(280), AB33 (360) 


colored flowers, 
late flowering 


I\/1QTL-Ae8/Ae9 [Ae-Ps3.1) 


III 


PI 180693 
552 


13 
10 


Field (RI-DI-LS-ATH), RRI-ADI 

Field (RI-DI-PLM), RRI-ADI CC 
(Rd84, Ael 09) 


8-27% 
6-14% 


22 - 54 


AA175 (280), AB92 (360), 
AD57 (320) 


late flowering 


MQTL-Ael2 {Ae-Ps4.1) 


IV 


552 


7 


Field (RI-DI), RRI-ADI 


8-21% 






- 






90-2131 


3 


Field (Rl), RRI-ADI 


8-13% 


6-42 


AD]86 (320), AA]74 (450) 








PI 180693 


3 


Field (RI-DI), RRI-ADI 


5-8% 








MQTL-Ael5 {Ae-Ps4.4, 
Ae-Ps4.5, Aphl) 


IV 


90-2079 


8 


Field (PLM-LS), RRI-ADI-DWL 
CC (SP7, Ael 06) 


6-44% 


8-28 


AC22 (200), AC32 (270) 


- 


MQTL-Ael6/Ael7 
{Ae-Ps5.1, Aph2) 


V 


90-2131 
90-2079 


7 
4 


Field (LS), ADI CC (RB84, Ae78, 
Ae85, Ae87 = SP7, Ael 06, Ael 09) 

Field (PLM, LS) CC (Ael 06) 


7-38% 
7-24% 


27-36 


AA81 (250), AB23 (390), AD68 
(290), AB47 (330) 


- 


MQTL-Ae25/Ae26 
[Ae-Ps7.6o) 


VII 


90-2131 


28 


Field (RI-DI-TPX-PLM-LS), 
RRI-ADI CC (RB84, Ae78, Ae85, 
Ael 06, Ael 09) 


6-60% 


39-54 


ADI 86 (320), AA174 (450) 
27-36 AA81 (250), AB23 (390), 
AD68 (290), AB47 (330) 39 - 54 


- 






PI 180693 


15 


Field (RI-DI-LS-ATH), RRI-ADI- 
DW CC (RB84) 


6-14% 




AA505 (180), AB27 (115), AB136 
(330), AB122 (330), AA387 (470), 
ABlOl (370) 








552 


10 


Field (RI-DI-PLM), RRI-ADI CC 
(Ael 09) 


6-20% 









^ Field locations/resistance evaluation criterion are coded as followed : Rl = Riec/Belon (FR); Dl = Dijon (FR), TPX = Templeux (FR), PLM = Pullman (USA); LS = LeSueur (USA); RRI = Root Rot Index; ADI = Aerial Decline 
Index; DW = Dwarfism; DWL=% Dried Weight Losses. Controlled-condition scoring traits are indicated by the name of the A. euteiches strain (RB84, Ae78, Ae85, Ae87 = SP7, Ael 06 and Ael 09). 
^ Range of percentage of phenotypic variation explained by each resistance-enhancing allele at each meta-QTL, depending on the variable. 

Confidence interval of the single or two meta-QTL region (1st value) and of the genomic region encompassing LOD-1 confidence intervals of all initial QTL underlying the meta-QTL (2nd value), in cM Kosambi. 
^ Recommended co-dominant SSR markers covering the genomic region encompassing LOD-1 confidence intervals of all initial QTL underlying the meta-QTL. Mean allele sizes are indicated in parentheses for 
each SSR. 



Hamon et al. BMC Plant Biology 201 3, 13:45 
http://www.bionnedcentral.conn/1471-2229/13/45 



Page 14 of 19 



resistance depending on the evolution of A. euteiches 
strains. Despite the limited potential for gene flow in soil- 
borne pathogens, A. euteiches could be assigned to the 
group of medium-moderately high evolutionary risk path- 
ogens (scale 5-7) [8], since the genetic structure of A, 
euteiches populations was reported to be patterned by a 
mixed reproduction system, including regular selfing with 
occasional migration of novel genotypes or outcrossing 
[74]. Consequently, pyramiding consistent Aphanomyces 
resistance QTL is strongly recommended, both for further 
increasing levels of resistance and for attempting to 
preserve the durability of each QTL effect. The best com- 
binations of QTL to increase resistance efficiency and 
durability have yet to be determined. Validation of QTL 
effects in different genetic backgrounds, knowledge of 
QTL effects on pathogen life cycle and epidemic develop- 
ment and an understanding of molecular mechanisms 
underlying QTL, will help identify QTL combinations that 
should be integrated in pyramiding strategies for durable 
resistance breeding. 

Methods 

Plant material 

A population of 111 Fio-derived pea recombinant inbred 
lines from the cross Dark Skin Perfection x 90-2131 was 
produced by single-seed descent at USDA-ARS, Pullman, 
WA, USA. DSP (Dark Skin Perfection, Unilever Limited 
1960, also designated as W6 17516) is a spring pea cultivar 
used for freezing and canning, with white flowers, normal 
leaves, wrinkled seeds with white hilum and it is suscep- 
tible to A, euteiches, 90-2131 [38] is a pea germplasm with 
white flowers, normal leaves, dimpled seeds with black 
hilum and it is partially resistant to A, euteiches in France 
and in the USA. The parental lines, 90-2131 and DSP, 
were used as check lines in field- and controlled-condition 
disease resistance tests of the DSP x 90-2131 RIL popula- 
tion. The three RIL parental lines previously studied in [46], 
Baccara (susceptible), PI180693 and 552 (partially 
resistant), were also used in the controlled condition 
assays. The MN313 line was included in the assays for 
distinguishing the two main pea pathotypes of A, 
euteiches [51]. The pea spring variety Solara (suscep- 
tible) was used as the adjacent control in the field 
assays of the RIL population. 

Pathogen material 

Six pure culture strains of A, euteiches virulent on pea 
were used for resistance evaluation of the DSP x 90- 
2131 RIL population in controlled conditions. These 
strains were RB84, Ael06, Ae78, Ae85, Ae87 (referred to 
as strain SP7 in [51]) and Ael09 (referred to as strain 
Ae467 in [75] and [33]). The six strains were chosen 
based on their various geographical origins and 



pathotype groups (Additional file 1). The strains Ael06 
and Ae87, were previously used for disease screening of 
the Puget X 90-2079 RIL population [45] and the strains 
RB84 and Ael09 were used in the Baccara x PI180693 
and Baccara x 552 RIL populations screenings [46]. 

Disease resistance evaluations 

In controlled conditions, resistance of the DSP x 90- 
2131 RILs to pure-culture strains of A, euteiches was 
assessed on 14-day-old seedlings in a climatic chamber, 
as described in [46], with one strain and four blocks per 
test. Each block included all RILs and control lines, with 
five plants/line/block. A root rot index (RRI) ranging 
from 0 to 5 was calculated as the mean disease severity 
score on the five plants per line, as described in [46] . 

In the field, the DSP x 90-2131 RILs were evaluated 
over an international Aphanomyces infested nursery net- 
work, described in [46], over four years (2000 to 2003) 
and five locations in France (Riec-sur-Belon, Finistere 
(RI); Dijon-Epoisses, Cote d'Or (DI); Templeux, Somme 
(TPX)) and in the United States (Pullman, WA (PLM); 
LeSueur, MN (LS)) [in 2000: PLM, LS; in 2001: LS; in 
2002 : RI, DI, LS; in 2003: RI, DI, TPX, PLM, LS]. Field 
assays were carried out using experimental designs as 
described in [46]. For field experiments conducted in 
France (two years, three locations) and in Pullman in 
2003, the design included a check plot of the susceptible 
cultivar Solara or DSP every two to four plots, in order 
to adjust disease severity scores for local disease varia- 
tions in the soil [76] using the formula described in 
[46,77]. Two disease criteria were used to assess resist- 
ance for each plot, as described in [46]: (i) the root 
rot index (RRI), using a 0-5 scoring scale, evaluated 
each year in French nurseries, except at DI in 2002, and 
(ii) the aerial decline index (ADI) evaluated once, twice 
or three times in all USA and French disease nurseries 
(except at RI in 2002), using a 1-5 and 1-9 disease scor- 
ing scale, respectively. 

Evaluation of morphological and phenological traits 

The DSP X 90-2131, Baccara x 552 and/or Baccara x 
PI 180693 RIL populations were evaluated for two agro- 
nomic traits: earliness at flowering (FLO) and plant height 
(HT). 

The FLO trait was evaluated in the Aphanomyces- 
infested nursery of DI in 2003 (DSP x 90-2131 RILs), 2004 
(Baccara x 552 RILs), 2006, 2007 and 2008 (Baccara x 
PI 180693 RILs) using the experimental design established 
for resistance evaluation. The FLO trait was also evaluated 
in a healthy nursery at Rennes-Le Rheu (lUe-et-Vilaine, FR 
(REN)) in 2002 (DSP x 90-2131 RILs), 2005 and 2008 
(Baccara x PI 180693 RILs) using a randomized complete 
block design with 1 and 3 block(s) in 2002/2005 and 2008, 
respectively (40 plants/plot in a two m-long twin rows). 
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The FLO trait was scored on each plot as the number of 
days to 50% bloom (FLOl) or to 100% bloom (FL02) from 
the first day of the year. The HT trait was evaluated in the 
same REN healthy nurseries as used for FLO evaluation in 
2002 (DSP X 90-2131 RILs), by measuring the average 
height of five plants at maturity in a whole plot. 

Molecular markers and genetic mapping 

The DSP X 90-2131 RIL population was genotyped using 
simple sequence repeat (SSR) from [50], random amplified 
polymorphic DNA (RAPD) [78] and known-function 
genes [48,49]. Two morphological traits were also scored: 
PI for hilum colour and R for round/wrinkled seeds. DNA 
extractions and PGR amplifications were performed as de- 
scribed in [46] and in [48,49,53]. The Puget x 90-2079 
RIL population was genotyped using additional SSR, 
RAPD, SCAR and known-function gene markers from 
[49,53,69,78-80] compared to markers reported in [31]. 
Marker coding and the genetic map (in cM Kosambi) 
from the DSP x 90-2131 RIL population were established 
with a minimum LOD score threshold of 3.0 and a max- 
imum recombination frequency of 0.4, as described in 
[46]. For each locus, adjustment of allelic segregation to 
the expected 1:1 Mendelian ratio was analyzed using a 
test (a=0.01). Additional markers were placed on the 
framework Puget x 90-2079 genetic map reported in [31], 
using the "assign" and "place" commands of MAP- 
MAKER/EXP version 3.0b [81]. 

Statistical and QTL analyses 

Statistical and QTL analyses were conducted from field 
and controlled-condition disease scoring data obtained in 
the DSP X 90-2131 RIL population and from earliness 
and plant height data obtained in the DSP x 90-2131, 
Baccara x PI 180693 and Baccara x 552 RIL populations, 
for each scoring variable in each environment. Statistical 
analysis of each data set was carried out using a two-way 
ANOVA estimating genotype and block effects and nor- 
mality of residuals was analyzed, as described in [46]. 
Broad sense heritability, RIL least-square means used for 
linkage analysis and Pearson correlation coefficients (r^) 
between adjusted mean data were also estimated as de- 
scribed in [46]. 

Additive-effect QTL analysis from each data set was 
performed by Composite Interval Mapping [82] using 
Windows QTL Cartographer 2.5 software [83], as de- 
scribed in [46]. Using the permutation test with 1000 
permutations, minimum LOD thresholds of 2.9 (for DSP x 
90-2131 RIL population) and 2.8 (for the Baccara x 
PI180693 and Baccara x 552 RIL populations) were 
chosen for all the traits to declare a putative QTL sig- 
nificant, corresponding to a genome-wide a error risk 
of 5%. Two QTL were considered as belonging to the 
same genomic region when their one-LOD drop-off 



confidence intervals overlapped. QTL for resistance to 
A. euteiches, earliness and plant height were named 
''Ae-Ps \ "Flo-Ps" and "HT'Ps \ respectively, followed by 
the linkage group number and the QTL number within 
the linkage group for each trait. Based on common 
markers between genetic maps, Aphanomyces resist- 
ance QTL common to those previously published 
[31,46] were named as described in [46]. The most sig- 
nificant pairwise epistatic interactions were detected for each 
resistance variable between all possible marker pairs of the 
DSP X 90-2131 genetic map, as described in [73] using a de- 
tection threshold of P < 7,1.10'^ and > 5%. 

From the Puget x 90-2079 RIL population, additive- 
effect QTL were re-detected, as described in [31], from 
the updated genetic map generated in this study and the 
phenotypic data reported in [31,45]. 

QTL meta-analysis 

Using all the additive -effect QTL identified for resistance 
to A. euteicheSy earliness and plant height from Puget x 
90-2079 [31,78], Baccara x PI180693, Baccara x 552 
[46] and DSP x 90-2131 (this study) RIL populations, a 
QTL meta-analysis was performed using the MetaQTL 
software version 1.0 [17]. 

QTL meta-analysis was conducted in three steps, 
according to details given in [21]. For each QTL, data 
given to the software were the QTL position (LG, position 
on the LG at the LOD peak), the upper and lower bound 
and LOD decrease of the QTL confidence interval, the 
percentage of variation [1^) individually explained by the 
QTL, the trait related to the QTL and the size of the map- 
ping population used for the QTL detection. First, a single 
consensus marker map was built by integrating the avail- 
able genetic maps, based on common markers designated 
with common names between the maps, using the 
ConsMap command of the software. The implemented 
method applied a weighted least-square strategy, using in- 
dividual distances between markers in each individual 
map, to determine marker order and position on the con- 
sensus marker map. The InfoMap command was carried 
out in order to list markers whose orders were not con- 
sistent between the different individual maps. Markers 
with inconsistent positions were removed from the con- 
sensus marker map. 

Second, the QTL detected in each study were projected 
onto the consensus map, using the QTLproj command of 
the software. This command enabled the homothetic pro- 
jection of the individual QTL positions and confidence 
intervals based on a scaling rule between QTL-flanking 
marker positions on the individual maps and on the con- 
sensus map. QTL projection was carried out using LOD-1 
confidence intervals of all individual QTL, for graphical 
representation and for identifying main genomic regions 
comprised of overlapping QTL intervals. For the meta- 
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analysis, the projection was carried out using independent 
individual QTL, since the QTL meta-analysis algorithm 
implemented in the software assumes that the input map- 
ping studies are independent from each other. Independ- 
ent QTL were selected as follows. In a given mapping 
population and for a given variable, when QTL detected 
from various environments or strains had overlapping 
confidence intervals, only the QTL with the greatest pro- 
portion of the phenotypic variation {R^) and the smallest 
confidence interval (if highest were equals) was retained. 
For example, among the five individual QTL detected at 
the top of LGII for the ADI variable in different environ- 
ments (RI04, RI05, PLM04, TPX04) in the Baccara x 552 
population, the one having the greatest R^ (15%) was kept 
for the meta-analysis and the four others were removed. 
Confidence intervals (CI) of independent QTL were esti- 
mated from the QTL R^ values, using the empirical formula 

proposed by [55]: CI = where N is the population size. 
This formula usually gives larger confidence intervals than 
the usual interval length of LOD-1 decrease. 

Finally, the QTL meta-analysis algorithm, computed in 
the QTLClust command of the software, was used to de- 
termine the most likely number of meta-QTL on a given 
chromosome and to estimate their corresponding posi- 
tions and confidence intervals. Meta-QTL considers QTL 
positions and corresponding confidence intervals in indi- 
vidual experiments, after projection onto the consensus 
map. We used the -cimode option 4 of QTLClust com- 
mand that consider the maximum confidence interval be- 
tween the LOD-1 decrease value reported and the QTL 
R^ value estimated by the formula of [55]. On a given 
chromosome, projected independent QTL were clustered 
into all possible numbers of hypothetic clusters of meta- 
QTL (/<), for which a Gaussian mixture model estimates 
meta-QTL positions and confidence intervals [17]. The 
optimal K was determined using the highest weight of evi- 
dence values estimated for five information-based criteria 
[17], which were computed for each K, The most frequent 
optimal /<r value given by the different criteria was selected, 
which always corresponded to the optimal K value deter- 
mined by the Akaike Information Criterion. The position, 
95% confidence interval and probability of individual QTL 
belonging to the meta-QTL, were given by the software 
for each meta-QTL. 

Two meta-analysis runs were conducted separately for 
Aphanomyces resistance QTL data and earliness and 
plant height QTL data, respectively, resulting in the de- 
tection of meta-QTL for two types of traits, which were 
then compared. 

Meta-QTL for resistance to A, euteiches and meta-QTL 
for earliness and plant height were named ''MQTL-Ae' 
and "MQTL-Morpho'\ respectively, followed by the meta- 
QTL number on the whole pea genome. 



Identification of candidate genes underlying meta-QTL 

Positional candidate genes included within the support 
interval of main resistance meta-QTL regions were identi- 
fied using the Pea-M. truncatula toolkit of Bordat et al 
[47]. Mining the high colinearity between the pea and M. 
truncatula genomes, the toolkit allowed placing in silico 
5460 pea Unigenes on the pea consensus map of [47], from 
positions of their best homologs on the M. truncatula gen- 
ome. Positions of main meta-QTL regions identified in this 
study were estimated on the pea consensus map of [47] 
based on common markers, mainly SSRs. On the pea map 
of [47], intervals between adjacent marker genes covering 
main meta-QTL regions were identified and a list of 
Unigenes contained in these intervals was established. In- 
tervals defined from adjacent marker genes genetically 
mapped in [47] were often larger than the reduced confi- 
dence interval of meta-QTL or absent. In these cases, inter- 
vals at most probable meta-QTL positions were chosen. 
Positional candidate genes were examined to identify those 
known to be involved in disease resistance in plants. 

Additional files 



Additional file 1: Strains oi Aphanomyces euteiches used, with code 
number, geographical origin and pathotype group. The pathotype 
groups of A euteiches strains were described ir^ Wicker ar^d Rouxel (2001) 
and defined by their differential reactions on a six pea genotypes (Wicl<er 
et al. 2003). 

Additional file 2: Frequency distribution of adjusted means of root 
rot index (RRI) scores for resistance to six strains of A.euteiches, in 
the DSP X 90-2131 pea RIL population. Values of the partially resistant 
(90-2131) and susceptible (DSP) parents, named Pr and Ps, respectively, 
are shown by arrows, '^n = total number of RILs assessed; ^m = mean ± 
standard deviation of the RIL population; ^h^ = heritability estimate. 
Strains: (a) RB84 (b) Ael06 (c) Ae85 (d) Ae87 (e) Ae78 (f) Ael09. 

Additional file 3: Frequency distributions of adjusted mean scores 
obtained in the DSP x 90-2131 RIL population for two 
Aphanomyces root rot resistance criteria (root rot and aerial decline 
indexes) assessed in 1 1 environments over four years and five 
locations in France and the USA. In each environment, one to three 
ADI scores were obtained. Adjusted mean values of the partially resistant 
(90-2131) and susceptible (DSP) parents, named Pr and Ps, respectively, 
are shown by arrows, = total number of RILs assessed; ^m = mean ± 
standard deviation of the RIL population; ^h^ = mean-based heritability of 
the trait. 

Additional file 4: Pearson correlation coefficients between the 
different adjusted mean scoring data obtained for Aphanomyces 
root rot resistance in the DSP x 90-2131 RIL population in 11 field 
environments and in controlled conditions with RB84, Ae106, Ae85, 
Ae87, Ae78 and Ae109 strains of A euteiches. Significance level 
threshold (P < 0.001) are indicated by asterisks {***). Field scoring 
variables are coded with abbreviations of the locations, years and scoring 
criteria tested: Rl = Riec/Belon (FR); Dl = Dijon-Epoisses (FR); TPX = 
Templeux (FR); PLM = Pullman (USA); LS = LeSueur (USA); 00 = 2000; 01 = 
2001; 02 = 2002; 03 = 2003; RRI = Root Rot Index; ADI = Aerial Decline 
Index. 

Additional file 5: Additive-effect QTL identified from the DSP x 90- 
2131 RIL population for resistance to A. euteiches in infected fields 
over 1 1 environments in France and in the USA, and in controlled 
conditions using six pure-culture strains of A euteiches. Values of 
QTL parameters were obtained using Windows QTL Cartographer 2.5 
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(LOD > 2.9). The QTL are ordered by positior^ or^ the LG. ^ Based or^ 
common markers, QTL common to Hamon et al. (201 1) were named as 
they were previously. ^ Field scoring traits are coded as follows: location 
(Rl = Riec/Belon (FR); Dl = Dijon (FR), TPX = Templeux (FR), PLM = Pullman 
(USA); LS = LeSueur (USA)); year (2000 = 00; 2001 = 01 ; 2002 = 02; 2003 = 
03); criterion (RRI = Root Rot Index; ADI = Aerial Decline Index). 
Controlled-condition scoring traits are indicated by the name of the 
strain (RB84, Ae78, Ae85, Ae87, Ael06 and Ael09). QTL position from 
the first marker of the linkage group (in centimorgans Kosambi) ^ Nearest 
left marker from the LOD score peak of the QTL ^ Log of likelihood ratio 
(LOD) value at the LOD peak of the QTL for each variable ^ Position of 
the lower and upper of the one-LOD drop-off confidence interval of the 
QTL, from the first marker of the linkage group (in centimorgans 
Kosambi) ^ Percentage of phenotypic variance explained by an individual 
QTL ^ Effect of substituting DSP alleles for 90 2131 alleles at the QTL. A 
positive sign indicates that QTL alleles increasing the resistance are 
contributed by the resistant parent 90-2131, whereas a negative sign 
means that resistant alleles are brought by the susceptible parent DSP. 

Additional file 6: Pairwise epistatic interactions associated with 
resistance to Aphanomyces root rot, identified from the DSP x 90- 
2131 RIL populations (P < 7,1.10'^). QTL are ordered by QTL name and 
decreasing R^ values. ^ The name of additive-effect QTL involved in 
epistatic interactions are indicated in bold. ^ Significance level of each 
pairwise epistatic interaction in an ANOVA model with two factors 
(markers) and an interaction factor (interaction between two markers) ^ 
Percentage of phenotypic variation explained by each individual 
interaction ^ Adjusted mean scores for the four genotypic classes defined 
by each marker pair : (SS) and (RR) DSP and 90-2131 alleles at the two 
markers, respectively; (SR) DSP allele at the first marker, 90-2131 allele at 
the second marker and conversely (RS). In bold: the most resistant 
genotypic class. 

Additional file 7: Frequency distributions of adjusted mean scores 
obtained in the DSP x 90-2131, Baccara x 552 and Baccara x 
P1 180693 RIL populations for plant height (HT, in cm) and flowering 
time (Flol and Flo2, in number of days from the first day of the 
year). Adjusted mean values of the partially resistant (90-2131, PI180693 
or 552) and susceptible (DSP or Baccara) parents, named Pr and Ps, 
respectively, are shown by arrows. °n = total number of RILs assessed; 
^m = mean ± standard deviation of the RIL population; = mean-based 
heritability calculated when scores were obtained on three blocks. 

Additional file 8: Additive-effect QTL identified from the DSP x 90- 
2131, Baccara x PI1 80693 and Baccara x 552 RIL populations for 
flowering time and plant height traits, scored in eight infested or 
healthy environments (values obtained by Windows QTL 
Cartographer 2.5, LOD > 2.9 for the DSP x 90-2131 RIL population 
and LOD > 2.8 for the Baccara x PI1 80693 and Baccara x 552 RIL 
populations). The QTL are ordered by position on the LG. Scoring 
traits are coded as follows: population (D2131 = DSP x 90 2131; B552 = 
Baccara x 552; BPI = Baccara x PI 180693); location (Dl = Dijon (FR); REN = 
Rennes (FR)); year (02 = 2002; 03 = 2003; 04 = 2004; 05 = 2005; 06 = 2006; 
01 = 2007; 08 = 2008); criterion (Flol and Flo2 = Flowering times; HT = 
Plant Height) ^ QTL position from the first marker of the linkage group 
(in centimorgans Kosambi) ^ Nearest left marker from the LOD score peak 
of the QTL Log of likelihood ratio (LOD) value at the LOD peak of the 
QTL for each variable ^ Position of the lower and upper of the one-LOD 
drop-off confidence interval of the QTL, from the first marker of the 
linkage group (in centimorgans Kosambi) ^ Percentage of phenotypic 
variance explained by each individual QTL ^ Effect of substituting Baccara 
or DSP alleles for 552 or Pll 80693 or 90-2131 alleles at the QTL. A 
positive sign indicates that QTL alleles increasing the earliness at 
flowering or decreasing plant height are contributed by the resistant 
parent 552, or Pll 80693 or 90-2131, whereas a negative sign means that 
earliness or short plant alleles are brought by the susceptible parent 
Baccara or DSP. 

Additional file 9: Localization of individual QTL and meta-QTL for 
resistance and morphological traits onto the consensus marker 
map: linkage group II, III, IV, V, VI and VII (for legend, see Figure 2). 

Additional file 10: Description of the pea consensus marker map 
established using the three individual genetic maps (Baccara x 



PI180693/552, DSP x 90-2131 and Puget x 90-2079) using the 
MetaQTL software. Cumulated genetic distances between markers are 
indicated in centimorgans Kosambi. The occurrence number of each 
marker across the three individual maps is indicated (Meta occurrence). 

Additional file 11: List of putative candidate genes underlying 
main Aphanomyces resistance meta-QTL regions, identified using 
the pea-M. truncatula translational toolkit [47]. Candidate genes and 
their contig sequence names were identified from positions of adjacent 
gene markers mapped on the consensus genetic map of [47], which 
were estimated to be located in the same regions as the meta-QTL by 
comparative mapping. 
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